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Abstract. We present Swift and XMM -Newton observations of the bright gamma-ray burst GRB 050326, detected by the Swift 
Burst Alert Telescope. The Swift X-Ray Telescope (XRT) and XMM-Newton discovered and the X-ray afterglow beginning 
54 min and 8.5 hr after the burst, respectively. The prompt GRB 050326 fluence was (7.7 ± 0.9) x 10~ 6 erg cnT 2 (20-150 keV), 
and its spectrum was hard, with a power law photon index T = 1.25 ± 0.03. The X-ray afterglow was quite bright, with a 
flux of 7 x 10~ u erg cnT 2 s -1 (0.3-8 keV), 1 hr after the burst. Its light curve did not show any break nor flares between 
~ 1 hr and ~ 6 d after the burst, and decayed with a slope a = 1.70 ± 0.05. The afterglow spectrum is well fitted by a power- 
law model, suffering absorption both in the Milky Way and in the host galaxy. The rest-frame Hydrogen column density is 
significant, Vh, z & 4x 10 21 cm -2 , and the redshift of the absorber was constrained to be z > 1.5. There was good agreement 
between the spatial, temporal, and spectral parameters as derived by Swift-X RT and XMM-Newton. By comparing the prompt 
and afterglow fluxes, we found that an early break probably occurred before the beginning of the XRT observation, similarly 
to many other cases observed by Swift. However, the properties of the GRB 050326 afterglow are well described by a spherical 
fireball expanding in a uniform external medium, so a further steepening is expected at later times. The lack of such a break 
allowed us to constrain the jet half-opening angle !?j ^7°. Using the redshift constraints provided by the X-ray analysis, we 
also estimated that the beaming-corrected gamma-ray energy was larger than 3 x 10 51 erg, at the high end of GRB energies. 

Despite the brightness in X rays, only deep limits could be placed by Swift- UVOT at optical and ultraviolet wavelengths. Thus, 
this GRB was a “truly dark” event, with the optical-to-X-ray spectrum violating the synchrotron limit. The optical and X-ray 
observations are therefore consistent either with an absorbed event or with a high-redshift one. To obey the Ghirlanda relation, 
a moderate/large redshift z S 4.5 is required. 
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1. Introduction 


The Swift satellite iGehrels et all 2004 ) is a mission dedi¬ 
e-mail: cated to the study of gamma-ray bursts (GRBs) and their af¬ 
terglows. GRBs are detected and localized by the Burst Alert 
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Telescope (BAT; Barthelmv et allEfXul) . and followed up atX- 
ray (0.2-10 keV) and optical/ultraviole t (1700-6000 A) wave- 
lengths by the X-Ray Telescope (XRT; iBurrows et alj 2005a I 
and th e Ultraviolet/Optical Telescope fUVOT: |Romina et al. 
12005al) . During the first year of operation. Swift has observed 
some 75 GRB afterglows, already doubling the pr e-Swift sam¬ 
ple. This rich dataset has allowed to study in detail the X-ray 
light curves, both at early and late timesGeadin<i to the dis 
covery of a complex behaviour (e.g iTagliaferrLeLaL L 2005a; 


covery ot a compiex benaviour (e.g agliarem et al., 2 
Nousek et all l2005t Chincarini et all 20051 ICusumano 


let al 


2005 3). Coupled with optical data, either from UVOT 
(e. g^Blusd n et all 20051 ) or ground-based observatories (e.g. 


Ber L>er*eTaiT 2005 ), this has opened a new era in the after¬ 


glow modeling. Swift also provided the first detection of truly 
dark GRBs, that is, events with no optical emission up to 
very deep limits dRoming et all l2005hl) . The study of high- 


burst at z > 6 (Watson et all 2005ti 

Cusumano et al. 

2005b 

Haislio et al. 2005; Price et all 200§ 

Tauliafcrri et al. 

2005b 

Kawai et all 2005). Moreover, it was found that Swift GRBs 


have a larger_average redshift than those discovered by earlier 
missions Ijakobsson et al 1 2005). 

During the performance verification and calibration phase 
(2004 Nov 20 through 2005 Apr 5), Swift observed sixteen 
GRB afterglows. Twelve of them were observed in automatic 
mode, and, among these, eight could be promptly (within 200 s 
since the trigger) observed by XRT and UVOT. In the remain¬ 
ing four cases, the beginning of the observation was delayed by 
approximately 50 min due to the Earth occultation constraints. 
This is the case for the bright GRB 050326, which was discov¬ 
ered by BAT on 2005 Mar 26 at 9:53:55 UT dMarkwardt et all 


2005). Its coordinates were arjiooo = 00 h 27 m 34 s , bnooo = 
71°22'34", with an uncertainty radius of 3' (95% contain¬ 


ment; Cummi ngs et all [2005 . Th is burst was also de tected by 
the Wind -Konus experiment Golenteskii et ilJ- • 2005), leading 
to the characterization of its broad-band gamma-ray spectrum. 

The Swift narrow field instruments could begin observing 
only 54 min after the BAT trigger. A bright, uncatalogued X-ray 
source was detected by XRT inside the BAT error circle, and 
was proposed to be the X-ray afterglow dMoretti et all 120053) . 
However, no source was detected by UVOT at this location 
^Holland et. all ’ 2005 ). XRT collected data up to 6.15 d after 
the burst. Subsequently, the decay of the light curve prevented 
any further detection of the afterglow. This object was_also ob¬ 
served for 45.8 ks b y XMM -Newton (Eh le & Perez Martinez! 
2005: De Luca el al.ll2005i . starting 8.5 hr after the trigger. 

Only limited ground-based follow-up was reported for this 
burst. This was likely due to its unfavorable location in the sky 
(very few telescopes can point at such low declination), as well 
as to the brightness of the Moon (which was 99% full at the 
time of the GRB explosion). No counterpart at wavelengths 
other than the X rays was reported. 

In this work, we present a complete discussion of the Swift 
and XMM-Newton observations of GRB 050326. In Sect.^Jwe 
describe the properties of the prompt emission. In Sect. 0 we 
describe in detail the XRT observations, the data reduction pro¬ 
cedure, and the temporal and spectral analysis; in Sect. 0 we do 
the same for the XMM-Newton data. In Sect.[5]we compare the 



Fig. 1 . The background-subtracted BAT light curve in the 20- 
150 keV energy band. The origin of the time axis was set to the 
instrument trigger, but a weak peak is apparent »9s before. 


results of the two instruments. In Sect.[6]we describe the UVOT 
optical observations. Finally, in Sect.^we present the physical 
implications of our observations in the framework of the stan¬ 
dard GRB afterglow model. Our conclusions are summarized 
in Sect. [8] 

Throughout this paper, all errors will be quoted at 90% con¬ 
fidence level for one parameter of interest, unless otherwise 
specified. The reduced^ 2 will be denoted as^ 2 , and the number 
of degrees of freedom with the abbreviation “d.o.f.”. We follow 
the convention F v (v, t ) oc t~ a v where a and [5 are the tem¬ 
poral decay slope and the spectral index, re spect ively. As time 
origin, we will adopt the BAT trigger (Mar kwardt et all 12005). 
The photon index is T = 1 + fi. Last, we adopt the standard 
“concordance” cosmology paramet ers, namely O m = 0.27, 
Q a = 0.73, h 0 = 0.71 (e.g. lSnergel et allEofU ). 

2. Prompt emission 

We reduced the BAT data using the latest available release of 
the HEADAS software (version 1.4). The light curve in the 
BAT energy band (20-150 keV) presents an initial, weak peak, 
9 s before the trigger, followed by several bright distinct peaks 
(Fig. 0. The Too and T$o durations of the burst (that is, the 
time intervals in which 90% and 50% of the fluence were col¬ 
lected, respectively) were 29.5 and 19.3 s, in the 20-150 keV 
band, respectively. We first modeled the BAT spectrum as a 
single power law with photon index T. This provided a good 
fit (x\ = 1.06 for 53 d.o.f.), yielding T = 1.25 ± 0.03 in the 
20-150 keV energy range. The fluence in the same band was 
(7.7 ± 0.9) x 10~ 6 erg cm -2 . 

GRB 05 0326 also triggered the W)W-Konus detector 
dGolenteskii et alJT 2005): in the 20 keV-3 MeV energy range it 
lasted 38 s, and had a flu ence of (3.22 + 0.05) x 10~ 5 erg cm~ 2 . 
iGolenteskii et al. ( 200 5) fitted the time-integrated spectrum of 
the burst as me asured by the Wind -Konus detector with a Band 
model (iBand et all 119931) . that is a smoothly joined broken 
power law with low- and high-energy photon indices Ti and 
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r 2 , respectively, and break energy Eq. The best fit provided 
Ti = 0.74 ± 0.09, T 2 = 2.49 ± 0.16, and E 0 = 160 ± 22 keV. 
The corresponding observed peak energy (that is, the energy at 
which the maximum of the emission is reached) was E p o \, s = 
(2 -Fy )£() = 200 + 30 keV. Motivated by their results, we also 
performed a fit to the BAT data using the Band function. Since 
the break energy Eq lies close to the upper boundary of the 
BAT energy range (150 keV), we were forced to freeze Eq and 
T 2 to the values determined by WiW-Konus. The fit was again 
good (xl = 1.14 for 53 d.o.f.), and provided Ti = 0.87 + 0.03, 
in go od agreement with the value found by iGolenteskii et alJ 
(20051. It is not surprising that both functional forms provide a 
good fit to the data, since they do not differ significantly inside 
the BAT energy range. Nevertheless, thanks to the very broad 
band covered by the VVznB-Konus instrument, for this burst the 
break energy could be clearly constrained. In the following, 
we will consider the Band model as the best description of the 
GRB 050326 spectrum. With this fit, the fluence in the 20-150 
keV band was (7.6 ± 0.8) x 10 6 erg cm~ 2 . Integrating the burst 
spectrum from 1 to 10 000 keV, we could evaluate the bolomet- 
ric fluence T of the burst, finding r F - 2.4 x 1 0 s erg cm~ 2 . 

No spectral evolution could be detected in the BAT data. 
We splitted the observation in three time intervals, covering the 
ranges [-9,-1], [-1,13], and [13,29] s (relative to the BAT 
trigger). By fitting the data with either a simple power law 
or with the Band model, the resulting parameters were always 
consistent with those derived by fitting the whole spectrum. 


3. XRT data analysis and results 

3.1. Data reduction 
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Fig. 2. XRT image of the field of GRB 050326, smoothed with 
a Gaussian kernel (477 full width half maximum). Events were 
accumulated from the first segment of the observation (589 s 
exposure time). The BAT refined position (black dot) is also 
shown together with its 95% containment error circle (3' ra¬ 
dius, dashed line). The XRT and XMM -Newton positions (177 
apart) are also plotted (crosses), and are almost indistinguish¬ 
able. The X-ray error circles (376 and 175 radius, respectively) 
are too small to be seen with this scale. 
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For a technical description of XRT an d its operations, we re¬ 
fer, e.g., to iBurrows et al ] ( 12005ah and (Hi 11 et al . (2004). XRT 
started observing the field of GRB 050326 on 2005 Mar 26 at 
10:48:27 UT, that is, 3307 s after the BAT trigger. The last ob¬ 
servation ended on 2005 Apr 1 at 13:30:53, i.e. 6.15 d after the 
burst. Occasionally, some reflected light from the Earth limb 
made the very low energy (< 0.2 keV) background increase 
abnormously, so that XRT incorrectly switched from the pho¬ 
ton counting (PC) to the windowed timing (WT) mode, even 
if the target count rate was well below 1 count s -1 . The ef¬ 
fective exposure time was 59.6 ks in PC mode and 20.7 ks 
in WT mode, leading to the collection of 614 and 580 pho¬ 
tons, respectively (0.3-10 keV energy band). As the satellite 
settled on the target, XRT recorded a source count rate of 
1.3 count s _1 , which dropped to 3 x 10~ 4 count s _1 at the end 
of the observing campaign (2005 Apr 1). From the third orbit 
after the start of the observation onwards, the source count rate 
was <0.1 count s -1 , while the background level was typically 
> 3 count s 1 over the whole field of view. Since WT data have 
only one dimensional spatial information, their S/N ratio was 
much lower than that of PC mode data. We therefore decided to 
consider WT data only for the first two orbits, when the source 
S/N was higher. Data were reduced using the xrtpipeline 
task of the latest available release of the HEADAS software 
(version 1.4). Accumulating the PC data from all observations, 
we found that the centroid position of the afterglow had coor¬ 


dinates a J2 ooo = 00 h 27 m 49716, fcooo = -71°22'1476, with a 
376 uncertainty radius (95% containment. Fig. 0. This posi¬ 
tion takes into account the correction for the misalignment be¬ 
tween the telescope and the satellite optical axis dMoretti et all 
l2005bl). This position is 1(3 a way from the refined BAT po¬ 
sition dCummings et al. , l2005b. a nd 374 aw ay from the pre¬ 
liminary XRT position i Moretti et alll2005al) . calculated using 
only the data from the first orbit and without the misalignment 
correction. 

3.2. Temporal analysis 

In order to extract the light curve, we considered all PC data, 
but discarded the WT data taken after the second orbit of the 
XRT observation (f > 20 ks). PC events were selected hav¬ 
ing grades 0-12 from a circle with 20 pixel radius (47"), cor¬ 
responding to 9 2% of the encircled energy fraction (EEF) at 
1.5 keV (Mor etti et alll2005cl) . Only the data in the 0.3-10 keV 
band energy range were considered (even if there are no events 
above 7 keV). To take into account the pile-up effect, during 
the initial part of the first orbit (f ^ 4000 s) an annular extrac¬ 
tion region with inner radius of 3 pixels (7") was adopted for 
PC data. This area includes 40% of the EEF, and the deriving 
PSF losses were consequently taken into account. The accu¬ 
racy of the PSF model in its central part is ~ 5% dMoretti et all 
l2005dl) . This error was properly propagated when evaluating 
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Time since trigger (s) + 10 s 


Fig. 3. The light curve of GRB 050326 and of its afterglow in the 0.3-8 keV energy band (see text for the computation of the flux 
conversion factors). XRT (black circles) and XMM -Newton data (empty diamonds) show a very good agreement (see also the 
inset). The solid line shows the fit to the combined XRT/XMM afterglow light curve. The dotted lines indicate the 90% errors of 
the extrapolated X-ray light curve. Light filled circles indicate the extrapolation of the BAT data to the 0.3-8 keV energy range, 
assuming the Band model as the best-fit spectrum. In this figure, the time origin was set 10 s before the nominal trigger time, 
to show the weak, untriggered precursor. This has no effect on the determination of the afterglow decay slope, due to the late 
beginning of the XRT observation. 


the final uncertainty of the PSF-corrected points in the light 
curve. The background in PC mode was evaluated by integrat¬ 
ing the signal from an annulus with inner and outer radii of 
50 and 90 pixels, respectively, centered at the afterglow posi¬ 
tion. Inside this region, the contamination from the afterglow 
is expected to be negligible. In WT mode, events were selected 
having grades 0-2 from a 20 pixel (47") wide rectangular re¬ 
gion, centered on the detector X coordinate of the afterglow. 
To estimate the background in WT mode, we considered a re¬ 
gion of the same size centered 40 pixels (94") away from the 
center of the afterglow. 

The XRT observation was split in different time segments 
because of the Earth occultation constraints. Each satellite orbit 
lasts ~ 5 800 s, while the target could be typically observed 
for approximately 1 000 s per orbit. To extract the light curve, 
the source events were binned in 10 s intervals, and these bins 
were further grouped to ensure a minimum of 50 counts per 
bin. When the counts in the last bin of each orbit were less 
than half of the required minimum (25 counts), the bin was 
merged with the previous one. From the fifth orbit onwards, 
XRT did not collect enough photons within a single orbit, so 
data from different orbits were merged. We eventually obtained 


a background-subtracted light curve composed by 25 points, 
with a minimum of 10 and a maximum of 57 counts per bin. 

The resulting light curve is shown in Fig. 0 and displays 
a uniform decay rate, with no indications of breaks or flares. 
A single power law fit provides a good description to the data, 
yielding a decay slope a = 1.64 + 0.07. In order to look for 
spectral variations across the observation, we computed the af¬ 
terglow hardness ratio as a function of time. To this end, we 
selected the events with energy below and above several piv¬ 
otal energies, and computed their number ratio. No significant 
variation was found during the whole observation, after setting 
the pivotal energy to 1, 1.5 and 2 keV. 

3.3. Spectral analysis 

To extract the spectrum of the source in WT mode, we used 
the same extraction regions, the same background regions, and 
the same screening as for the temporal analysis. For PC mode, 
however, we further selected only grade CM- events in order to 
improve the spectral resolution. The spectrum was binned in 
order to ensure a minimum of 20 counts per energy bin, ig¬ 
noring channels below 0.3 keV. The spectral analysis was per¬ 
formed using XSPEC (vll.3). We first considered WT and PC 
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Fig. 4. The spectra of the GRB 050326 afterglow as observed by XMM-Newton (PN: squares; MOSl: open circles; MOS2: 
triangles) and by Swift-X RT (filled circles). The solid lines represent the best-fit absorbed power-law model convolved with the 
instrumental responses (see Table 13 for the best-fit parameters). XRT data were selected from the time interval covered by the 
XMM-Newton observation. 


data separately, using the former for the first 2 orbits (3 to 15 ks 
after the burst) and the latter for the rest of the observation. In 
both cases, the spectrum was fitted with an absorbed power- 
law model, yielding good^f 2 values (Table [3 . The best-fit val¬ 
ues for the Hydrogen column density A+ and for the photon 
index T did not show significant variations between the first 
(WT data) and second part (PC data) of the observation. In 
fact, combining the data from the two segments together, we 
obtained an excellent^ 2 value, indicating that the spectral prop¬ 
erties of the afterglow did not change during the observation. 
Kg. a shows the XRT spectrum (filled circles), together with 
the best-fit absorbed power-law model (precisely, only data si¬ 
multaneous to the XMM-Newton observation were used for the 
plot; see Sect.|3. 

To look for the presence of absorbing material in 
the proximity of the afterglow, we tried to estimate the 
Galactic Hydrogen column density /Vh.mw towards the 
GR B direction. We fou nd three different measurements: 
iDickev & Lockamnl(ll99()l) give (Vh.mw = 4.6 x IQ 20 cm 2 ; the 
Leide n/Argentine/Bonn Galactic HI Survey JKalherla et all 
120051) provides a l ower value, JVhmw = 3.8 x 10 20 cm -2 ; last, 
the dust maps by dSchlegel et all 11993) give A v = 0.12 mag, 
which corresponds to (Vhm w = 2.2 x IQ 20 cm~ 2 after assum¬ 
ing the prescription given bv IPredehl & Schmitt ( 1995 1. While 
the average of these three independent measurements is (3.5 + 
1.7) x 10 20 cm~ 2 , we conservative lyadoptcd the largest of the 
above values (4.6 x 10 20 cm~ 2 l. lStratta et al.l ( 2004 ) estimated 


Table 1. Best-fit spectral parameters for the two segments of 
the XRT observation (WT and PC data), fitted both separately 
and together. We report the results either leaving the Hydrogen 
column density as a free parameter, or freezing it to the Galactic 
value. 



N h (10 21 cm- 2 ) 

r 

xi Cr 2 /d.o.f.) 

XRT (WT) 

1 3 +0 - 7 

1 ‘ -0.6 

1 q 9+0.27 

1 -0.24 

1. 10(29.1/27) 

XRT (PC) 

1 4+O.8 

1 '^-0.7 

2 00 + °- 38 
z,uu -0.31 

1.05 (12.6/12) 

XRT (WT+PC) 

1 4+O.6 

1 -0.5 

1 95 +0 - 21 

1.10(42.0/38) 

XRT (WT) 

0.46 (frozen) 

1.58 + 0.12 

1.22 (34.1/28) 

XRT (PC) 

0.46 (frozen) 

1.58 + 0.14 

1.41 (19.1/13) 

XRT (WT+PC) 

0.46 (frozen) 

1.61 ±0.10 

1.20(47.0/39) 


that th e typical error affecting the maps by Dic kev & Tockamnl 
dl99(j is 30%, which is not far from the scatter among the three 
measurements. The best-fit Hydrogen column density derived 
by the XRT afterglow spectrum is marginally unconsistent with 
the Galactic value. Fixing /V H = /Vh.mw provided a poor fit 
(xl = 1-20 for 39 d.o.f.). The probability of such a worsening 
in the fit is < 7.5%, as estimated by an F-test. Therefore, XRT 
data marginally suggest the presence of additional absorbing 
material, likely located in the GRB rest frame. In the next sec¬ 
tion, we will present further evidence for the presence of excess 
absorption, based on XMM-Newton data with better S/N ratio. 
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4. XMM-Newton data analysis and results 

The afterglow of GRB 050326 was also observed by XMM- 
Newton as a target of opportunity, starting on 2005 Mar 
26 at 18:25 UT (8.5 hr after the burst). The observation 
lasted for 45.8 ks. Data were collected with the European 
Photon Ima ging Camera (EPIC), which consists of the 
PN dStriider et all 12001 1) and of the two MOS detectors 
iTurner et all 1200 li All the cameras were operated in full- 
frame mode with a thin and medium optical filter on PN and 
MOS, respectively. A preli minary analysis of these data was 
presented bv lDe Luca et alJ ( 2005a l. 


4.1. Data and temporal analysis 

The appropriate observation data files were retrieved from 
the XMM Science Archive. The data reduction was per¬ 
formed using the most recent release of the XMM Science 
Analysis Software (SAS v6.1.0), with the standard pipeline 
tasks (epprocand emprocforPN and MOS, respectively). The 
observation was badly affected by high particle background 
(soft proton flares), with almost no nominal (quiescent) back¬ 
ground time intervals. The back-illuminated PN CCD is partic¬ 
ularly sensitive to this background; indeed, more than 25% of 
the PN observing time was lost due to the detector switching 
to its counting mode 1 . The afterglow of GRB 050326 (source 
XMMU J002748.8-712217 :iEhle & Perez Martinez! EPOS ) was 
anyway clearly detected in all cameras. The astrometry of the 
EPIC images was improved by cross-correlating serendipitous 
X-ray sources in the field with objects in the USNO-B1 cata¬ 
log. This yielded the following refined coordinates for the af¬ 
terglow: Qrj2ooo = 00 h 27 m 4971, c5 J2 ooo = -71°22'1673, with an 
1-cr uncertainty of 175. The EPIC and XRT positions differ by 
177, and are therefore fully consistent within the uncertainties. 

In order to retain a S/N ratio large enough to perform the 
temporal and spectral analysis, a standard time-filtering ap¬ 
proach to screen soft proton flares could not be applied (a high 
particle flux was present during the whole observation). Thus, 
source events were extracted with a particularly stringent spa¬ 
tial selection, considering only the innermost portion of the 
point spread function. We used a circle of 15" radius (con¬ 
taining ~ 65% of the EEF). The PSF correction was applied 
to the flux and spectral measurements by computing the ad-hoc 
effective area using the SAS task arf gen. The error in this pro¬ 
cedure is estimated to be at most at the 5% level 2 and it was 
properly taken into account in the light curve error budget. 

Background events were selected from source-free regions 
within the same CCD chip where the source was imaged. In 
particular, for the PN data we used 2 boxes of 45" x25" located 
at the same distance from the readout node as the target; for the 
MOS we used an annulus centered at the target position with 
inner and outer radii of 90" and 180", respectively. With such 


1 The counting mode is activated when the count rate in a quad¬ 
rant exceeds the telemetry limit (~ 400 count s _1 for the PN). In this 
mode, the informations for individual events of that quadrant are not 
transmitted to ground. 

2 See page 9 of 

http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0818-2- 


a choice, the background amounted to ~ 13% and ~ 9% of the 
counts in the source extraction region for PN and MOS data, 
respectively, in the 0.3-8 keV range. The overall (background- 
subtracted) number of source events was 3 990,1 850 and 1 760 
in the PN, MOS 1 and MOS2 detectors, respectively. 

The background-subtracted count rate clearly showed a de¬ 
clining trend with time. We again fitted the light curve assum¬ 
ing a power law decay. The value of the decay slope a was eval¬ 
uated independently using PN, MOS 1 and MOS2 data, yield¬ 
ing fully consistent results. We therefore repeated the fit us¬ 
ing the combined dataset, finding a = 1.72 ± 0.09 in the 0.3- 
8 keV energy range (xl = 1.30 for 40 d.o.f.). The background- 
subtracted light curve is shown in Fig-EJ together with the XRT 
light curve, after converting the count rates to fluxes using the 
best-fit absorbed power-law models described in Sect. I3.3l and 

IO >. 

4.2. Spectral analysis 

Spectra for the source and the background were extracted from 
the same regions used for the temporal analysis, as described 
above. Source spectra were rebinned in order to have at least 
30 counts per energy bin and to oversample the instrumental 
energy resolution by a factor 3. Ad-hoc response matrices and 
effective area files were created with the SAS tasks rmfgen 
and arfgen, respectively. The spectral analysis was performed 
using XSPEC (vl 1.3). The spectra were fitted simultaneously 
in the 0.3-8 keV band. Since the MOS observation started ~ 
1 hr earlier than the PN, a PN/MOS normalization factor was 
introduced in the fit as a further, free parameter. Due to the 
fading of the source, this also implies that the observed time- 
averaged flux is expected to be higher in the MOS than in the 
PN. 

An absorbed power-law model reproduced the spectrum 
quite well (xl = 1.17 for 190 d.o.f.). The best-fit parameter 
values are reported in Table[2](lirst row); the MOS-PN normal¬ 
ization factor was 1.08+0.04. Fig.@]shows the spectra collected 
by the EPIC cameras together with the best-fit model. Both the 
photon index and the Hydrogen column density are in good 
agreement with those found by XRT, but are much better con¬ 
strained. In particular, the value of Ah inferred from the fit was 
significantly larger than the Galactic one Ah.mw (Sect. ED. 
Moreover fixing Ah = Ah.mw resulted in a much poorer fit 
(xl = 1-85 for 191 d.o.f.), and even increasing Ah,mw by 30% 
(see Sect. 13.31 and Istratta et all I2OO4I ) the fit was still unac¬ 
ceptable (xl = 1.65 for 191 d.o.f.). A significant improvement 
was achieved by fixing Ah = Ah.mw and adding to the spec¬ 
tral model an extra neutral absorber at redshift z with column 
density Ah, z . This yielded^ = 1.03 (189 d.o.f.); the chance 
probability of such improvement, as estimated by an F-test, was 
< 1 x 10 6 with respect to the model containing only one ab¬ 
sorbing component (with free Ah). After adding the extra ab¬ 
sorption component, the best-fit power law photon index was 
T = 2.03 + 0.05, while the intrinsic gas column density and the 
redshift were X H .- ~ 6 x 10 22 cm -2 and z ~ 6, respectively. 
However, the latter two values are not well constrained, ow- 
4 .ipifto their strong correlation (Fig.H}. In any case, the spectral 
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Fig. 5. Confidence contours (68%, 90% and 99% levels for 2 
parameters of interest) for the gas column density A H .- and the 
redshift z of the intrinsic absorber, as computed from the fit to 
the EPIC spectra. The Galactic colum n density was assum ed to 
be Ah.mw = 4.6 x 10 2 ° cm 2 (iDickev & Lockamnlll990l) . The 
inset shows a zoom-in of the low-redshift region. 

fit allowed us to constrain Ah,; > 4 x 10 21 cm~ 2 and z > 1.5 
(at 90% confidence level for 2 parameters of interest; see inset 
in Fig. |3- We investigated the dependence of these confidence 
contours on the assumed value of the Galactic column density. 
By varying Ah.mw by 50% (within the range discussed above), 
we found that the 90% confidence interval on Ah,- varies by 
20%, while that on z varies by less than 10%. 

The observed (time-averaged) fluxes in the 0.2-10 keV 
band were ~ 5.1 x 10~ 13 and ~ 5.5 x 10~ 13 erg cm~ 2 s -1 in 
the PN and MOS, respectively. The corresponding unabsorbed 
fluxes were ~ 7.4 x 10~ 13 and ~ 7.9 x 10~ 13 erg cm~ 2 s _1 . 

In order to search for possible line features in the spectrum 
(both in emission and in absorption) we divided the 0.4-5 keV 
energy range into 0.2 keV intervals. For the continuum, we as¬ 
sumed the best-fit absorbed power-law model (including also 
the rest-frame column density). For each of the intervals, we 
added a Gaussian line with fixed width (smaller than the in¬ 
strumental resolution) and central energy free to vary within 
the selected interval; the normalization could be either posi¬ 
tive or negative. We then repeated the exercise with different 
choices of the energy intervals. We found no significant lines 
in the 0.4-5 keV range in the combined MOS/PN dataset. The 
upper limit (3-cr) on the equivalent width of any line is ~ 50 
and ~ 250 eV in the 0.4-2 and 2-5 keV energy ranges, respec¬ 
tively. 

We also tried to fit the spectr a with thermal models (e.g. 
[Reeves et all 2002 ; £ azzat 112003). A redshifted, optically thin 
plasma emission model (MEKAF in XSPEC) was used, with the 
redshift linked to that of the intrinsic absorber. The fit worsened 
ixl ~ 1.4) with respect to the simple power law model, either 
fixing the metal abundances to Solar values, or leaving them as 
free parameters. 

Finally, we looked for possible spectral evolution with time, 
using the power law plus redshifted absorber model described 


Table 2. The spectral parameters as measured separately by 
XMM -Newton and XRT, fitting the data both separately and 
together. For this fit, XRT data were selected from the same 
time interval covered by the XMM -Newton observation. Due 
to the limited statistics, Ah was frozen to the value derived by 
XMM -Newton when fitting XRT data. The errors are at 90% 
confidence level for a single parameter of interest. 


Instrument 

A h (10 21 cm- 2 ) 

r 

x 2 y (r/d.o.f.) 

EPIC 

XRT 

EPIC + XRT 

1 3+ 0 ' 1 

1.3 (frozen) 

1 3 +01 

o i q+0.06 

^ -0.06 
i 06+0.29 

1 -0.27 

? 09 +0 05 
^•^-0.08 

1.17 (222.5/190) 
1.43 (11.4/8) 

1.18 (219.7/197) 


above. For this study we divided the data into two subsets with 
exposure times of ~ 15.8 and ~ 27.6 ks, each subset containing 
approximately half of the afterglow counts. We then extracted 
the corresponding spectra for the source and the background. 
No significant (> 3-cr) variations in the spectral parameters 
were found (except, of course, for the flux normalization). 

5. XMM-NewtonISwift-XRJ comparison 

As noted above, the spectral parameters as derived using XRT 
and XMM -Newton data nicely agree within the errors. To bet¬ 
ter check the consistency between the afterglow temporal and 
spectral properties as derived by the two satellites, we per¬ 
formed a more accurate operation. We selected XRT data from 
the MOS observing time interval (35-76 ks after the burst). 
During this period, the XMM -Newton observation was contin¬ 
uous, whereas Swift completed eight orbits, providing 13 ks of 
effective exposure time for XRT. Therefore, in order to com¬ 
pare the two datasets, we had to assume that both the spectrum 
and the light curve behaved in a regular fashion. For example, a 
flare during a Swift occultation would distort the results of the 
comparison. This assumption seems anyway fully justified by 
the XMM -Newton data, which show a regular behaviour of the 
light curve without significant irregularities. 

We fitted XRT data using the absorbed power law model. 
For the sake of comparison between the two instruments, we 
did not add any rest-frame absorption component. Due to the 
limited statistics, moreover, we froze the column density to the 
value found by XMM-Newton (which is fully consistent with 
that found by XRT for the full observation). The results of the 
fit are presented in Table |3 Within the errors, the photon in¬ 
dex is consistent with that found by XMM -Newton. Last, we 
fitted together the XRT and XMM -Newton data from the three 
EPIC instruments (leaving Ah as a free parameter), obtaining 
our final best fit (Table[2}. 

Our second step was to compare the flux normalization fac¬ 
tors for the two instruments. To this extent, we selected the 
data only from exactly overlapping intervals (8.5 ks effective 
observation time). We froze Ah and T to the values found pre¬ 
viously (which rest on a better statistics), computing only the 
normalization factors (the ratio of the fluxes observed by Swift- 
XRT and the EPIC detectors). The fit provided 0.90 + 0.14 for 
MOS1, 0.95 + 0.16 for MOS2, and 0.89 ± 0.15 for PN. This 
result indicate that, within the errors, the XRT current absolute 
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flux calibration is good, providing perhaps slightly underesti¬ 
mated values (at the ~ 5% level). 

As shown in the previous sections, the light curves from 
XRT and XMM -Newton were well described by power laws 
with slopes 1.64 + 0.07 and 1.72 + 0.09, respectively. These val¬ 
ues are consistent within their errors. Moreover, selecting XRT 
data from the 8.5 ks time interval with simultaneous XMM- 
Newton observations, we found a slope a = 1.73 ± 0.12, in 
perfect agreement with the XMM -Newton result (and consis¬ 
tent with that measured by XRT for the whole observation). 

In order to directly compare the fluxes measured from the 
two satellites, we converted the 0.3-8 keV count rates to unab¬ 
sorbed fluxes in the same band (Fig. 0. To compute the con¬ 
version factor for the XMM -Newton spectrum, we used the ab¬ 
sorbed (one component) power law model with the parameters 
reported in Table |3(first row). For XRT, we calculated a con¬ 
version factor for each of the two operational modes, using in 
both cases the absorbed power law best-fit model (leaving Ah 
as a free parameter), as reported in Tabled (first and second 
rows). The combined fit of the joint XRT and XMM data pro¬ 
vided a — 1.70 ± 0.05 (Fig. [3}. In the following, we will adopt 
this value as the best determination of the temporal decay slope. 


6. Optical and ultraviolet observations 


The UVOT instrument onboard Swift observed the field of 
GRB 050326 together with XRT, starting 54 min after the trig¬ 
ger. In the subsequent orbits, it collected a series of images in 
its 6 broad-band filters (V, B, U, UVW1, UVM2, and UVW2; 
Table QJ. The afterglow was not detected in any of the sin¬ 
gle or coadded exposures. Summing the images in each of the 
six filters, we estimated the 3-cr upper limits using the UVOT 
dedicated software (task uvotsource). The counts were ex¬ 
tracted from a 6" and 12" radius aperture for the optical and ul¬ 
traviolet filters, respectively, after subtracting the background. 
We then corre c ted th e upper limits for Galactic absorp tion 

f chlegel et all 49981 assuming the extinction curve of iPeil 
992l) . Our final limits are summarize d in Table |3l With re - 
spect to the original values reported by iHollandetal ( 2005 ). 
our measurements were obtained adopting the most recent on- 
flight calibration. 

The only reported ground-based optical observation for this 
burst was an R-band upper limit provided by the 0.6m telescope 
at the Mt. John Observatory (tlristram et al. i.; 2005 ). This mea¬ 
surement is also listed in Tabled 


iNousek et al. . 2005 ). lasting up to 10 3 —10 5 s (and in some 
cases even longer). Then, the light curve steepens again, lead¬ 
ing to a « 1-1.5; this phase was the one seen by BeppoSAX, 
XMM -Newton and Chandra. At late times,_a further steepen¬ 
ing is sometimes observed (e.g. VaughjmeUd,, 1200 5). likely 
the signature of a jetted outflow SE oa dsi ll999l) .In some cases, 
bumps and flares appear superimposed to the power-law decay, 
up to several tens of ks (e.g. Bur rows et alll2005bl) . 


The light curve of GRB 050326 exhibited a different be¬ 
haviour with respect to that outlined above. Its light curve 
showed a single, unbroken decay from » 55 min to » 4.2 d. 
However, our coverage began relatively late, so that we may 
have missed early deviations from the power law behaviour. 


In order to investigate the afterglow early stages and to an¬ 
alyze the connection between the prompt and afterglow emis¬ 
sion, we extrapolated the afterglow flux to the time of the 
prompt emission. We then compared the obtained value with 
that expected from the prompt emission in the XRT band (0.3— 
8 keV), computed adopting the Band best-fit model. The result 
is shown in Fig.0 where the light circles indicate the prompt 
emission fluxes. Since the GRB spectrum is known in good de¬ 
tail (particularly, since no breaks are expected between the BAT 
and XRT ranges), the extrapolation process should be quite re¬ 
liable. As it can be seen, if no temporal breaks were present 
in the X-ray light curve, the afterglow flux in the X-ray range 
exceeded the prompt one by a factor of ~ 100 (with a small un¬ 
certainty, due the tiny error in the decay index). We cannot ex¬ 
clude that such emission was present (since we have no prompt 
observations in the X-ray band), but, if present, the present 
component would appear as a very brigh t, soft excess . Such 
feature would not be unprecedented (Vandersn ek et all 12004: 
Ivetere et al 1120051) . but in this case it would likely contaminate 
the low-frequency end of the BAT spectrum. Moreover, the soft 
excesses always contained less energy that the GRB proper. 
Thus, the most conservative hypothesis is to assume that a 
break was present in the early light curve, or that the after¬ 
glow onset was delayed. Indeed, as mentioned above, most of 
Swift afterglows show a shallow decline phase during the first 
thousands seconds after the GRB. Independently of any extrap¬ 
olation, we note that GRB 050326 was distinctly different from 
most bursts observed by BeppoSAX, for which the backward 
extrapolation of the late-time X-ray afterglow roughly matched 
the prompt emission level in the X-ray range, as measured by 
the Wide Field Cameras iFrontera et allEfiofl) . 


7. Discussion 

7.1. The X-ray light curve 

To date. Swift has observed X-ray emission from dozens of 
GRB afterglows. A systematic analysis of their light curves 
has evidenced several common features JNousek et all 12005: 
IChincarini et all 120051) . During the first few hundred seconds, 
a steep decay is often observed (a ~ 3-5; Tagli aferri et all 
12005 al) . us ually interpreted as the ta il emission from the prompt 
GRB ('e.g. lCusumano et~ahl 12005 al) . T his phase is followed bv 
a much flatter decline (a ~ 0-0.7; e.g. lCamnana et all 12005 a: 


We also performed a different operation. Using the best- 
fit X-ray spectrum, we extrapolated the XRT flux to the BAT 
energy range (20-150 keV), and reported it at the time of the 
burst using the afterglow decay law. Also in this case, the ex¬ 
pected value exceeded the observed prompt emission, but by 
a smaller factor. This again suggests that a break in the light 
curve was present before the beginning of the XRT observa¬ 
tion, but the evidence is less compelling. For example, we can¬ 
not even exclude that the afterglow spectrum had a break be¬ 
tween the XRT and BAT ranges, so that the extrapolation actu¬ 
ally overestimated its flux. 
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Table 3. Summary of optical and ultraviolet observations. All data but the point in the R band were measured by Swift-FWOT. The 
bv lTristram et al. (. 2005 1. All measurements were corrected for the Galactic extinction (Ay = 0.123 mag; 
assuming the Milky Way extinction curve bv IPei ( 199 2). The optical-to-X-ray spectral index is fox = 
\og(F° vt /F*)/ log(vx/v op t)- The X-ray frequency was set to vx = 1.5 keV, the logarithmic mean of the XRT observing range. 
The X-ray flux was computed by interpolating the XRT light curve to the time of the optical limit. The last column reports the 
ratio of the extrapolated X-ray flux to the optical one (assuming no spectral breaks between the two bands; see Fig. |7J. 


/Gbaud observations i: 
Schlegel et aL 19981) . 


Time since burst 

Filter 

Wavelength 

Exposure time 

Extinction 

Magnitude 

Flux density 

ftox 

X-ray/optical ratio 

(s) 


(A) 

(s) 

(mag) 


(/By) 



3300 

UVW2 

1930 

10 

0.29 

> 17.41 

<99.11 

< 0.52 

> 

3306 

V 

5460 

100 

0.12 

> 18.66 

< 107.5 

< 0.45 

> 65 +3 j 

3350 

UVM2 

2220 

100 

0.36 

> 18.54 

< 34.17 

< 0.32 

> 75/2i 

3470 

UVW1 

2600 

100 

0.27 

> 18.75 

< 27.80 

<0.29 

> 103/« 

5219 

u 

3450 

750 

0.20 

> 19.65 

< 11.00 

< 0.24 

> 176 

9390 

B 

4350 

714 

0.16 

> 20.87 

< 17.47 

<0.46 


24840 

Rmoa 

6399 

600 

0.10 

> 20.20 

< 25.70 

< 0.74 

> io/j" 


7.2. Constraints on the afterglow parameters 


The properties of the explosion can be in ferred in the con- 

model (e.g. iMeszaros & Reesi 
context, the observed emission 
is due to synchrotron radiation from a decelerating relativis¬ 
tic shock, which produces a decaying flux with a power-law 
spectrum. Depending on the model parameters, definite rela¬ 
tions between the spectral and temporal indices a and ft are 
predicted. The combined XMM -Newton and Swift-XKT data 
provide a = 1.70 + 0.05 and ft - 1.09 + 0.08. Both values are 


text of the standard afterglow 
In this 



2005|. These numbers are consistent with a spherical outflow 


expanding inside a homogeneous medium, if the XRT range 
was between the injection and cooling frequencies (vi and v c , 
respectively). In this case, the model prediction is a - 3/1/2 = 
1.63 ± 0.12, in excellent agreement with the measured value 
a = 1.70 ± 0.05. All other possibilities (a wind-stratified 
medium, or a different location of the break frequencies) are 
excluded at > 3.5-cr level. The power-law index of the elec¬ 
tron energy distribution is p — 1 + 4o-/3 = 1 + 2ft , so that 
p = 3.25 ± 0.06. Such value is rather high, but not unprece¬ 
dented. 


No break was observed in the X-ray light curve of 
GRB 050326 between 55 min and ~ 4.2 d after the burst. The 
condition Vj < v < v c thus held during this time range. While 
Vi typically lies below ~ 10 15 Hz for t > 1 hr (e.g. I.Sari et all 
1998^ Panaitescu & Kumail 2000 1). keeping v < v c up to t > 


4.2 d requires noe ^-2 < 3 x 10 11 ft, where n = no x 1 cm 


-5 i 


- 1/2 


is the ambient particle density, eg = 10 2 Sb -2 is the mag¬ 
netic field energy fraction, and E IS0 = £i s0 ,54 x 10 54 erg is the 
(isotropic-equi valent) fireball energy. This condition is difficul t 
to satisfy (e.g. IPanaitescu & Kuman l200lt lYost et all 120031) . 
so it may be regarded as a problem for the model. We note, 
however, that for this burst both a and ft could be measured 
with good accuracy, so the consistency between the predicted 
and observed value of the decay index is remarkable. 


The absence of any break also poses some constraints on 
the geometry of the emission. GRB afterglow light curves of¬ 
ten show a late-time steepening, commonly interpreted as the 


result of a jetted geometry. For GRB 050326, this break likely 
occurred after the end of the Swift observations. In fact, breaks 
earlier than ~ 1 hr are usually due to different reasons (suc h 
as the end of the refreshed shock episode; e.g. IZhang et all 
12005 ). If interpreted as a jet break, a very narrow jet would 
be implied. Moreover, the decay slope in the monitored time 
range is quite flat compared with t hat expected (and usually ob¬ 


served) after jet breaks ( e.g. Ilsrael et al. 


pected (an 

lll999tlH 


arrison et all 


h00~f: Klose et all 120041) . The measured decay would imply a 
hard electron distribution (p < 2). Using the relations provided 
bv lDai & Cheng! (l200ll) . an unreasonably low p — 4a - 6 = 
0.8 ± 0.2 would result (with the X-ray band being above v c ), 
which would give rise to a spectrum completely inconsistent 
with the observed one. Assuming v c above the observed range 
would only worsen the situation. 

Assuming that the jet break occurred at to ^ 4 d allows 
us to put some constraints on the jet opening angle and on the 
beaming-corrected gamma-ray energy. To this extent, the main 
obstacle is the lack of a well constrained redshift. Our spectral 
analysis of the XMM -Newton data allowed only to set a broad 
range 1.5 < z < 8 (90% confidence level). Using the bolometric 
fluence F of the prompt emission, the lower limit on z provides 
a constraint on the GRB radiated energy: 


DHz) 


= 4n-±-—F > 1.4 x 10 53 erg. 


( 1 ) 


1 + z 

where D L > 3.40 x 10 28 cm is the luminosity dis¬ 

tance. GRB 050326 was therefore quite likely bright in 
gamma rays_ compared to other GRBs detected by Swift (e.g. 

IChincarini et al lEooi. r _, 

Following lSarietai . ( 1999b . the jet half-opening angle can 
be constrained as follows: 


> 6.7° 


fii £ r 3/8 fi^[ 1 V^) 1/8 i (2) 

\ 2.5 / \10 53 erg/ \0.21cm 3 / 


where r/ is the prompt radiative efficiency and we have assumed 
fb > 4 d. Note that the dependance from the gamma-ray energy 
and on the other parameters is rather mild. The corresponding 
beaming-corrected gamma-ray energy is then 
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Fig. 6. Comparison of GRB 050326 with the Amati (right) and Ghirlanda (left) relations (jAmatij 
120041 ). The thick solid curves (black and grey) show the position of GRB 050326 as its redshift varies in the interval 0.1 < z < 10. 
The Ghirlanda track is actually a boundary (as the horizontal arrows indicate), since we can infer only a lower limit to the 
beaming-corrected energy at each redshift. Filled circles and squares indicate the GRBs which define the above two relations, 
plotted as straight solid lines_(together with their J-, 2- and 3-crcontours: long-dashed, dot-dashed and dotted lines, respectively). 
Data were taken from iGhirlanda et ah d2004E00 5T Grey diamonds indicate the intersection of the GRB 050326 tracks with the 
3-cr contours of the Amati and Ghirlanda relations. These points thus define the 3-cr redshift ranges for which GRB 050326 was 
consistent with the two relations. In the two GRB 050326 tracks, the region 1.5 < z < 8 (indicated by the X-ray data) is shown in 
black, bound by asterisks. 


The quantity Z\/( 1 + z) has only a mild dependence upon 
Z. For the redshift range allowed by our X-ray measurements 
(1.5 < z < 8), we have 1.4 x 10 28 cm < D L /( 1 + z) < 
2.8 x 10 28 cm, so that E y j > (3-8) x 10 51 erg. This is at the 
high end of the beaming-corrected gamma-ray energy distribu¬ 
tion (iGhirlanda et al.ll2004l) . 

Fig. [6] shows the position of GRB 050326 in the plane E y 
vs £p, to chec k how it c ompares with the Amat i and Ghirlanda 
relations dAmati et all l2002t IGhirlanda et all 120041) . As the 
redshift varies, the burst position follows the tracks indicated 
by the thick solid curves (see the figure caption for the de¬ 
tails). The Ghirlanda track is actually a boundary, since the 
lower limit on translates into a lower limit for E y ,j for any 
given redshift. We can ask for which redshifts GRB 050326 
was consistent with the two above relations. Only loose limits 
are provided by the Amati relation (which has a large disper¬ 
sion): at the 3-cr level, z > 0.25 is implied. The comparison 


with the Ghirlanda relation is less solid, since further assump¬ 
tions are needed (such as the ambient particle density and the 
break time). However, for our fiducial values, we have two al¬ 
lowed ranges: a low-redshift region (0.1 ^ z $ 0.8), plus a 
high-redshift solution (z <: 4.5). A lower particle density would 
move the track towards the left, while a larger jet break time 
would shift it rightwards. 

We note that the only way to make GRB 050326 in agree¬ 
ment with the Ghirlanda relation and simultaneously satisfy our 
constraints on the X-ray absorption (z ^ 1.5) is to require a high 
redshift for this event, since the low-redshift region is excluded 
by the fit to the X-ray column density at > 99% confidence 
level. Therefore, although our arguments are rather specula¬ 
tive, and would surely need more conclusive data, we regard 
GRB 050326 a s a m oderate/high redshift candidate. 

Recently, Liang & Zhang (;2 005l ) have presented a model- 
independent multidimensional correlation between the ob- 
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served isotropic energy £ r i so , the rest-frame spectral peak en¬ 
ergy £p, and the comoving break time f h /(l + z). This relation, 
in principle, allows to compute E yps0 (and z.) for a GRB with 
known E p and f h . However, no significant constraints could be 
inferred in the case of GRB 050326. Moreover, since this re¬ 
lation was derived using the break time as measured in the 

to X-ray data may not be valid 


optical band, its application 
fLiang & Zhang 


7.3. Evidence for intrinsic absorption 


The presence of intrinsic absorption, besides allowing us to 
constrain the GRB redshift, has other important consequences. 
The rest-frame absorbing system has a Hydrogen column den¬ 
sity larger than ~ 4 x 10 21 cm~ 2 . For moderate redshifts, Nh.z 
would be much larger. Several afterglow observations, both 
from Swift and previous missions, showed evidence for ex¬ 
cess X-ray absorption in addition to t h e Gal a ctic value (e.g. 
Galama & Wiierst |200l|_ Stratta et^k[ 2004i iDe Luca et all 
2005 b|). Recently, ICamTOnaeTalTHoOSbl) showed that about 
half of Swift afterglows have a large rest-frame column den - 
sity, typical of giant molecular clouds dReichart & PriceH2002l) . 
Given the connection between GRB explosion and supernova e 
dGalama et all fl99i Istanek et a ll Eooi iHiorth rtlfl Eool. 


this fact may constitute a powerful way to study the regions 
where m assive star formation tak es place in the high-redshift 
Universe. lLazzati~& Pernald2002l) showed that the prompt GRB 
flux is able to ionize the surrounding medium up to radii 
as large as ~ 5 pc, therefore leaving no absorbing material. 
Such process may have been observed in act for GRB 000528 
dFrontera et all 12004 ). The fact that a large column density was 
measured in GRB 050326 may imply that the absorbing mate¬ 
rial was distributed in a wide region (R ^ 5 pc), or that the 
ionizing flux was not large. 


For GRB afterglows, comparable absorption in the X-ray 
range is usua lly not accompanied by large extinction in the 
optical band dGalama & Wiiersl l200ll) . Only small amounts 
of dust are usually inferred from the analysis of optical spec¬ 
tra, even when heavy extinction is observed in the X-ray af¬ 
terglow. Several explanations were invoked to explain this 
discrepance, among which the destructio n of dust from the 
burst and/or afterglow photons dWaxman & Drain tl 120008), a 
large gas-to-dust ratio in the intervening material (j^tratt^eUff , 
2Q24J), or an overabundance of a elements dWatson et al. . 
2005 a). Unluckily, several factors hamper the study of this 
problem, such as the uncertaintes in the shape and normaliza¬ 
tion of the extinction curve, the possibility that GRBs occur 
in special, low-metallicity e nvironments dFvnbo et all |2003; 
McFadven & Wooslev, 1 999j), and in several cases the lack of 
the redshift determination. 


For GRB 050326, no detection in the optical/ultraviolet 
band could be obtained. TablcQJrcports the available upper lim¬ 
its to the afterglow flux, from both Swift-UW OT and ground- 
based observations. In Fig. 0 we show the optical-to-X-ray 
spectral energy distribution at different epochs. We computed 
the X-ray flux at the time of each available limit, adopting 



Fig. 7. Broad-band spectral energy distribution of the afterglow 
of GRB 050326, computed at different times (which are iden¬ 
tified by different symbols). The shape of the X-ray spectrum 
was assumed to be constant throughout the observation, and the 
decay law was adopted to report the X-ray flux at the time of 
the optical measurements. 


the decay law measured by Swift-X RT and XMM -Newton. 
Furthermore, no spectral evolution was assumed. 

As it can be seen, the UVOT limits provide strong 
constraints, even if they are not particularly deep. In fact, 
GRB 050326 was bright in X rays (with a flux of ~ 
10 10 erg cm~ 2 s 1 , 1 hr after the GRB). Table 0 reports the 
optical-to-X-ray spectral indices pox ■ For all our measure¬ 
ments, Pox < Px- Moreover, most of them violate the syn¬ 
chrotron limit Pox > 0.5, which holds for a non-obscured 
synchrotron spectrum. GRB 050326 can therefore be classi¬ 
fied as truly dark, according to the definition proposed by 
[lakobsson et al ■ (2004). This limit is quite robust, since no as¬ 
sumptions are made about the position of the synchrotron break 
frequencies. For this GRB, however, we can go further in this 
reasoning. Our analysis of the temporal and spectral properties 
of the afterglow, in fact, indicated that the XRT range was be¬ 
low the cooling frequency v c . Therefore, the extrapolation of 
the XRT spectrum to the optical domain seems in this case reli¬ 
able, since no spectral breaks are expected to lie between these 
two bands. This allows estimating the suppression factor suf¬ 
fered by the optical flux, and is reported in the last column of 
Table 0 Again, large lower limits were found, implying con¬ 
spicuous rest-frame extinction (up to a factor ~ 100 and more, 
corresponding to > 5 mag). The presence of the injection fre¬ 
quency Vi close to or blueward of the optical band may partly 
expl ain the flux dearth. How ever, following the formulation of 
IPanaitescu & Kumar] d2000l) . even choosing rather extreme pa¬ 
rameters (£ is0 = 10 54 erg, z — 5, s e — sb — 0.1), v, can at most 
be comparable to the ultraviolet observed frequencies. In par- 
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ticular, at the time of the U -filter measurement, which provides 
the strongest constraint, v; cannot be blueward of this band. So, 
even if Vj has some role in this game, it cannot be responsible 
for the whole suppression of the optical flux. Moreover, low 
values of e B were required to keep the cooling frequency out¬ 
side the XRT range (see Sect. O. so that vj was likely at much 
lower energies than the optical band. 

The truly dark nature of this burst allows one of the follow¬ 
ing two possibilities. The burst may have suffered dust extinc¬ 
tion in its host galaxy. The amount of dust is not straightfor¬ 
ward to evaluate. The main obstacle is again the lack of the 
redshift, together with the unkwown shape of the extinction 
curve. However, our limit that more than 5 mag were missing in 
the observed U band may roughly correspond to A v ^ 2 mag 
for z ~ 1.5, even if many other solutions are acceptable. The 
second possibility is that GRB 050326 was at high redshift, as 
suggested by our analysis of the X-ray spectrum combined with 
the limits provided by the Ghirlanda relation. In this case, vir¬ 
tually no flux is left blueward of the redshifted Lyman dropout. 
To suppress the flux in the V band, z £ 5 would be required. 
However, the combination of a moderate redshift and mild ab¬ 
sorption may relax this condition. 

8. Conclusions 

We have presented a detailed analysis of the GRB 050326 
prompt and afterglow emission. The combined capabilities of 
Swift (which sampled the light curve for a relatively long time 
span) and XMM-Newton (which ensured a large statistics) al¬ 
lowed to obtain a thorough characterization of the afterglow 
properties. 

The prompt emission was relatively bright (with a 20- 
150 keV fluence of ~ 8 x 10 r ' erg cm -2 ). The spectrum was 
hard (photon index T = 1.25 + 0.03), suggesting a peak energy 
at the high end of the BAT energy range or beyond. Indeed, 
thanks to the simul taneous detection of this burst by the Wind- 
Konus experiment (iGolenteskii et al .:, 1200 5). the prompt spec¬ 
trum could be fully characterized. The prompt bolometric flu¬ 
ence was T ~ 2.4 x 10 -5 erg cm -2 (1-10000 keV), and the 
observed peak energy was £p i0 b s = 200 + 30 keV. 

Due to pointing constraints, XRT and UVOT observations 
could start only 54 min after the GRB. The X-ray afteglow was 
quite bright, with a flux of 7 x 10 -11 erg cm -2 s -1 (0.3-8 keV) 
1 hr after the GRB. However, no optical counterpart could be 
detected. The X-ray light curve showed a steady decline, with 
no breaks or flares. The best-fit power-law decay index was 
a = 1.70 + 0.05. Such regular behaviour is different from that 
usually observed by Swift, but this may be the result of the lim¬ 
ited time coverage (observations could be carried out only be¬ 
tween 54 min and 4.2 d after the burst). Indeed, extrapolation 
of the afterglow light curve to the time of the prompt emission 
overpredicts the burst flux, and may suggest a slower decay be¬ 
fore the beginning of the XRT observation. 

The analysis of the combined XRT and XMM-Newton data 
allowed to characterize in detail the afterglow spectrum. A fit 
with an absorbed power-law model provided a good description 
to the data, yielding a photon index T = 2.09 ± 0.08 and a col¬ 
umn density significantly in excess to the Galactic value. The 


best-fit model was thus computed adding an extra absorption 
component, leaving its redshift z free to vary. Although both 
Nu, z an d Z could not be effectively constrained, a firm lower 
limit Ah , z > 4xl0 21 cm -2 could be set. Therefore, GRB 050326 
adds to the gro wing set of afterglows with large rest-frame 
column density JGalama & WiiersL l200lt IStratta et all [2004; 
ICamnana et alll2005bl) . The limits measured in the optical and 
ultraviolet region by UVOT lie well below the extrapolation of 
the X-ray spectrum. In particular, they violate the synchrotron 
limit that the optical-to-X-ray spectral index should be larger 
than 0.5. This implies a large extinction and/or a high redshift. 

The X-ray spectral analysis also allowed us to set the lower 
limit z > 1.5 to the redshift of the absorbing component (and 
therefore of the GRB). The isotropic-equivalent gamma-ray en¬ 
ergy was then E y - K0 > 1.4 x 10 53 erg. The temporal and spectral 
properties of the afterglow were nicely consistent with a spher¬ 
ical fireball expanding in a uniform medium, with the cool¬ 
ing frequency above the X-ray range. We could therefore set a 
lower limit to the jet break time ~ 4 d. The jet opening angle 
could be constrained to be i) ] k. 7°, with only a weak depen¬ 
dence on the (unknown) fireball energy. The beaming-corrected 
gamma-ray energy was E y j = (3-8) x 10 51 (fb/4 d) 3 / 4 erg, in¬ 
dependently from the redshift. GRB 050326, thus, released a 
large amount of gamma_ray_s_(only GRB 990123 had a larger 
energy in the sample of iGhirlanda et all 12004 1. 

To be consistent with the Ghirlanda relation 
icihirlanda et all 200 4.1 two redshift ranges are allowed, 
either at low (z ^ 0.8) or high (z <; 4.5) redshift. However, to 
simultaneously satisfy also the limits derived from the X-ray 
spectral analysis, only the high-redshift region is left. We note 
that the Ghirlanda relation is still based upon a small sample, 
so that any inference cannot yet be regarded as conclusive. 
However, the results from the X-ray spectra, the consistency 
of the GRB 050326 properties with the Ghirlanda relation, 
and the strong dearth of optical/ultraviolet afterglow flux, are 
overall consistent with a moderate/high redshift (z ^ 4). A 
search for the host galaxy through deep infrared and optical 
imaging may conclusively settle this issue. 
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